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ABSTRACT 

Young solar-type stars rotate rapidly and many are magnetically active. Some appear to undergo magnetic 
cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating 
suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the 
convective envelope of a solar- type star rotating at 5 times the current solar rotation rate. We find that dynamo 
action builds substantial organized global-scale magnetic fields in the midst of the convection zone. Striking 
magnetic wreaths span the convection zone and coexist with the turbulent convection. A surprising feature of 
this wreath-building dynamo is its rich time-dependence. The dynamo exhibits cyclic activity and undergoes 
quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a 
roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow 
and are more organized temporally and spatially than those realized in our previous simulations of the solar 
dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our 
simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily 
through an 0-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily 
well represented by a simple a-effect. During a reversal the magnetic wreaths propagate towards the polar 
regions, and this appears to arise from a poleward propagating dynamo wave. As the magnetic fields wax and 
wane in strength and flip in polarity, the primary response in the convective flows involves the axisymmetric 
differential rotation which varies on similar time scales. Bands of relatively fast and slow fluid propagate 
towards the poles on time scales of roughly 500 days and are associated with the magnetic structures that 
propagate in the same fashion. In the Sun, similar patterns are observed in the poleward branch of the torsional 
oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. 
Subject headings: convection - MHD - stars:interiors - stars:rotation - stars: magnetic fields - Sun:interior 



1. INTRODUCTION 

Motivated by the rich observational landscape of stellar 
magnetism and by our own Sun's cycles of activity, we have 
undertaken 3D simulations of convection and dynamo action 
in solar type stars. These simulations have explored how 
global-scale flows of differential rotation and me ridional cir- 
culati on are established in rapidly rotating suns ( Brown et alj 
120081) and how dynamo action m ight be realized in these stars. 

In Paper I (B rown et al.l 1201 0) we studied dynamo action 
in a star rotating three times faster than our current Sun and 
found that global-scale magnetic structures can arise naturally 
in the midst of the stellar convection zone. These wreaths of 
magnetism were stable for long periods of time, maintaining 
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their identity for many thousands of days, and did not require 
a stably stratified tachocline to survive the turbulent convec- 
tion. Rather they coexisted with it. 

In this paper we explore a convective dynamo simulation 
in a solar-like star rotating five times the current solar rate. 
Here, global-scale magnetic wreaths still form in the convec- 
tion zone but now they become time dependent and undergo 
repeated cycles of magnetic polarity reversal. These wreaths 
are significantly more complex than the steady wreaths dis- 
cussed in Paper I. 

L L Observational Landscape 

Magnetism is a nearly ubiquitous feature of solar-type 
stars. Young, rapidly rotating suns appear to have much 
stronger magnetic fields at their surfaces. Observations re- 
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veal a clear correlation between rotation and magnetic activ- 
ity, as inf erred from proxies such as X-ray and chro mospheric 
emission (I Schriiver & Zwaarj2000l:|Pizzolato et al.i2003i) . but 
the observational landscape is complex with few ot her well- 
established trends to constrain dynamo models (e.g.. lRempell 
I2OO8; Lanza 2010). 

A fundamental characteristic of the solar dynamo is its reg- 
ular cycles of sunspots. Indeed, the 22-year solar activity 
cycle stands out as one of the most remarkable and enig- 
matic examples of magnetic self-organization in nature. So- 
lar type stars appear to undergo magnetic cycles with periods 
ranging from several years to several decades. Most obser- 
vational investigations of stellar activity cycles rely on long- 
term monitoring of photospheric, chromospheric, or coronal 
emission using the Sun as a baseline for calibration and com- 
parison (e.g., Baliunas et al. 1995; Hempelmann et al. 1996; 
Messina & Guinan 2002; Hall 2008; Olah et al. 2009; Lanzal 
2010^ Photometric and spectroscopic variability are linked to 
magnetic activity and periodic modulation on time scales of 
years to decades is interpreted as an activity cycle. The most 
extensive such survey is still the Mount Wilson HK Project 
which monitored chromospheric emission in 111 solar-type 
stars from 1966-1991 (BaHunas et al. 1995). About half of 
the stars in the sample (51) show clear signs of cyclic activity, 
including 21 with well-defined cycle periods ranging between 
7 and 25 years. The remaining stars exhibit either irregular 
variability with no clear systematic variation (29) or smooth 
time series with flat or linear trends (31). Longer-term but 
more sporadic photometric measurements are available for a 
few dozen stars and these, along with the Mount Wilson data, 
show evidence for multiple periodicities and possible varia- 
tions in the apparent starspot cycle pe riods over the course 
of multiple decades dOlah et al.ll2009t) . More rapidly rotat- 
ing stars generally have shorter cycles (ISaar & BrandenburgI 
[1999), but robust scaling relationships with stellar mass, ro- 
tation rate and other fundamental parameters have remained 
elusive. The shortest stellar activity cycles observed to date 
have periods of rou ghly 1 to 1 .6 years, and occur in rapidly 
rotating F-type stars (iMetcalfe et al.l20TotlGarcia et al.l2()10l) . 

1 .2. Elements of Cyclic Activity 

The magnetic fields observed at the surfaces of late-type 
stars must arise from dynamo action in the stellar convection 
zones, as in the solar dynamo. Self-organization in turbulent 
dynamos is intimately associated with helicity and shear In 
many mean-field dynamo theories, poloidal magnetic fields 
are built by turbulent correlations in the convection through 
what is known as an a-effect, and these correlations are en- 
hanced if the flow is helical. Rotation and stratification impart 
helicity, both kinetic and magnetic, and can also lead to large- 
scale shearing flows of differential rotation. Global-scale 
shear promotes self-organization through a process known as 
the fi-effect where mean toroidal fields are generated from 
poloidal fields. Dynamos employing these two effects are 
known as a-Vl dynamos, and these models form the basis of 
much of our theoretical understanding of cyclic dynamos is 
stars like the sun (e.g., Charbonneau 2010). 

The central role of shear and helicity in establishing cyclic 
activity has been confirmed by global numerical simula- 
tions of convective dynamos. The first 3D magnetohydrody- 
namic (MHD) simulations of convection in rotating spherical 
shells to exhib it cyclic behavior were explored in detail by 
iGilmanI (119831) . Recent simulations of global convective dy- 
namos have focused on rapidly-rotating Boussinesq systems 



in deep convective shells, often in th e context of planetary dy- 
namos (reviewed by^ usse 2000; Ro berts & Glatzmaieri200lt 
IChristensen & Aubertll2006 ). Here the magnetic field is often 
dominated by the axisymmetric dipole component and is typ- 
ically either stable in time or undergoes chaotic reversals. In 
planetary dynamos, rapid rotation, deep shell geometries, and 
minimal density stratifications promote quasi-2D convective 
columns that are strongly aligned with the rotation axis. Stel- 
lar convection occurs under quite different conditions. 

In the solar convection zone, density stratification is of fun- 
damental importance and the rotational influence is moder- 
ate (Rossby number Ro ^ 1), giving rise to intricate, three- 
dimensional, highly turbulent conv ective structures spa nning 
many spatial and temporal scales dMiesch et al.ll2008h . Dy- 
namo simulations in this parameter regime produce complex 
magnetic topologies, with more than 95% of the magnetic en- 
ergy in the fluctuating (non-axisymmetric) field components 
(iBrun et al.' 2004). Mean magnetic fields are likewise com- 
plex, with multipolar structure and transient toroidal ribbons 
and sheets. Polarity reversals of the dipole component occur 
but they are irregular in time. 

The presence of an overshoot region below the convec- 
tion zone and a tachocline of rotational shear there pro- 
motes mean-field generation, producing persistent bands of 
toroidal flux antisymmetric about the equator while strength- 
ening and stabihzing the dipol e moment jBrowning et al.l 
120061 120071: iMiesch et al. l '2009^. However, these simula- 
tions did not exhibit sys tematic magnetic cycles. Recent 
results by (Ghizaru et al. 2010) achieve both large-scale or- 
ganization and cyclic reversals of polarity in a convection 
simulation with an underlying tachocline. In that simula- 
tion, substantial mean magnetic fields are present in both the 
tachocline and the bulk of the convection zone. The simu- 
lation spans roughly 93,000 days (225 years) with polarity 
reversals occurring on roughly 11,000 day (30 year) periods. 
Convective simulations in spherical wedge geometries with 
self-consistently established differential rotation profiles and 
imposed tachoclines also achieve organi zed fields in the con - 
vection zone and cyclic dynamo activity jKapvla et alj|20ldl) . 
Interestingly, large-scale magnetic fields and organized polar- 
ity reversals have also been achieved in spherical geometries 
with helically forced flows, and involving neither glob al-scale 
differential rotation nor tachoclines (Mit ra et al .1120 10b . 

Here we explore convection and dynamo action in a rapidly 
rotating sun which spins five times faster than the current so- 
lar rate, calling this case D5. This dynamo builds globally 
organized fields but also undergoes quasi-regular reversals of 
magnetic polarity. This is achieved in the stellar convection 
zone itself; a tachocline is not included in our simulation and 
thus can play no role. We describe briefly how our simulations 
are conducted in Section |2] and then compare cyclic case D5 
with the dynamo case D3 from Paper I in Section|3] In Section 
IHwe examine the nature of convection in these rapidly rotat- 
ing suns. The magnetic fields and their global-scale polarity 
reversals are explored in Sections |5]-|8] The differential rota- 
tion shows marked signatures of the reversals, and these tor- 
sional oscillations are discussed in Section |9] Unusual mag- 
netic structures arise during certain intervals, and these are 
briefly discussed in Section [TO] We reflect on our findings in 
Sectionfm 
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2. SIMULATING STELLAR CONVECTIVE DYNAMOS 

We study MHD stellar convection and dynamo action with 
the anelastic spherical harmonic (ASH) code ( Clune et"aD 
119991: iBrun et al.ll2004l) . Our simulation approach is briefly 
described here, but the equations, implementatio n and nature 
of ou r study are explained in detail in Paper I ( Brown et al.l 
l2010h . This global-scale code simulates a stratified spher- 
ical shell and here we focus on the bulk of the convection 
zone, with our computational domain extending from Q.IIRq 
to 0.97 /?Q (with Rq the solar radius) and spanning about 
172 Mm in radius. We avoid regions near the stellar surface 
and also do not include the stably stratified region near the 
base of the convection known as the tachocline. The total den- 
sity contrast across the shell is about 25, with a stratification 
derived from a ID solar structure model (IBrown et al.ll2008l) . 
As discussed in Paper I (Section 3.1), we chose to simplify 
these calculations by omitting a tachocline of penetration and 
rotational shear, motivated also by the lack of observational 
evidence for tachoclines in rapidly rotating stars. We have ex- 
plored some companion simulations that include a tachocline, 
and find that the wreath-building dynamics and their temporal 
behavior are largely unmodified by its presence. We will re- 
port on this in another paper Furthermore, recent mean-field 
models of the solar dynamo suggest that the latitudinal shear 
in the lower convection zone is primarily responsible for the 
generation of the mean toroidal field, as opposed t o the ra- 
dial shear i n the tacho cline (Dikpat i & Gilmanll2006l: iRempell 
12006; .Munoz-Jaramillo et al. 2009) . The role of a tachocline 
in establishing cyclic magnetic activity remains unclear. 

ASH is a large-eddy simulation (LES) code with subgrid- 
scale (SGS) treatments for turbulent diffusivities. As in 
Paper I the vorticity, entropy and magnetic field diffusivities, 
V, K and rj respectively, are functions of radius only and vary 
with the background density as /T^/^. Thus, all three diffusiv- 
ities are smallest in the lower convection zone. As in Paper I, 
the magnetic Prandtl number Pm =v/ri = 0.5. The fundamen- 
tal characteristics of our simulations and parameter definitions 
are presented in Table [1] 

The dynamo simulation (case D5) was initiated from a ma- 
ture hydrodynamic progenitor which had been evolved for 
more than 8000 days and was well equilibrated (case H5). 
To initiate our dynamo case, a small seed dipole magnetic 
field was introduced and evolved via the induction equation. 
The energy in the magnetic fields is initially many orders of 
magnitude smaller than the energy contained in the convec- 
tive motions, but these fields are amplified by shear and grow 
to become comparable in energy to the convective motions. 
Our magnetic boundary conditions are a perfect conductor at 
the bottom of the convection zone and match onto an external 
potential field at the top. 

Stellar dynamo simulations are computationally intensive, 
requiring both high resolutions to correctly represent the ve- 
locity fields and long time evolution to capture the equili- 
brated dynamo behavior, which may include cyclic variations 
on time scales of several years. The strong magnetic fields 
can produce rapidly moving Alfven waves which seriously 
restrict the Courant-Friedrichs-Lewy (CFL) timestep limits in 
the upper portions of the convection zone. Case D5, rotating 
five times faster than the current Sun, has been evolved for 
over 17000 days (or over 8 million timesteps). As a historical 
note, with the higher resolution of this simulation, this rep- 
resents roughly a factor of a m illion more computation than 
was possible in lGilmanl ( ll983i) . which is in surprisingly good 
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Figure 1. Magnetic wreaths achieved in (a) case D3 (Paper I) and (h) case D5 
(this paper). Shown are time-latitude plots of mean toroidal (longitudinal) 
magnetic field (B^) at mid-convection zone, with scaling values indicated. 
Case D3 builds persistent, time-independent wreaths but the wreaths achieved 
in case D5 undergo quasi-regular reversals of polarity (three shown here, with 
roughly a 1500 day period). The dynamic range of the color bars are indi- 
cated. 

agreement with Moore's law doubling over the almost thirty 
year interval separating these simulations. We plan to report 
on a variety of other dynamo cases, some at higher turbulence 
levels and rotation rates, in subsequent papers; cyclic activity 
is a generic feature of many of these dynamos. 

3. DYNAMOS IN RAPIDLY ROTATING SUNS 

In order to accentuate and investigate the self-organization 
processes associated with helicity and rotational shear, we 
have conducted a series of simulations of solar-type stars ro- 
tating more rapidly than the Sun. The non-magnetic ana- 
logues of our dynamo simulations exhibit a systematic in- 
crease in the rotational shear of differential rotation with in- 
creasing rotation rate. Convection remains vigorous at all 
rotation rates, though at the highest rotation rates novel lo- 
calized nests of convection arise in the equatorial regions 
(Brown et al. 2008). 

In Paper I (Br own et al.l2 010) we describe in detail the gen- 
eration of persistent toroidal wreathes in a simulation rotating 
at three times the solar rate (case D3, at 3 ri©). These wreaths 
are localized bands of strong (^ 7 kG, with peak amplitudes 
of ^ 26 kG) toroidal flux located in the midst of a turbulent 
convection zone, sustained by rotational shear. The toroidal 
field strength peaks near the base of the convection zone at 
latitudes of ±15°, with opposite polarity in the northern and 
southern hemispheres. The mean toroidal (longitudinal) mag- 
netic fields for case D3 are shown in Figure [!};. These mag- 
netic structures form within 2000 days from weak initial seed 
fields and persist for the remainder of the simulation, span- 
ning over 15,000 days of evolution. They are nearly steady 
in time and do not show global-scale reversals of magnetic 
polarity. 
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Table 1 

Parameters for Primary Simulations 



Case 



Ra 



Ta 



Re Re' Rm Rm' 



Ro 



Ro' 



Roc 



D5 
H5 



96x256 x512 
96x256x512 



1.05x10* 
1.27x10* 



6.70x10' 
6.70x10' 



273 
576 



133 
141 



136 66 



0.273 
0.303 



0.173 
0.182 



0.241 
0.268 



0.940 
0.940 



1.8 



D3 
H3 



96x256x512 
96x256x512 



3.22x10' 1.22x10' 
4.10x10^ 1.22x10' 



173 
335 



105 86 
105 — 



52 0.378 
— 0.427 



0.255 
0.265 



0.311 
0.353 



1.32 
1.32 



2.64 



Note. — Dynamo simulation at 5 times the solar rotation rate is case D5, and its hydrodynamic (non-magnetic) companion is H5. 
Both simulations have inner radius rbot = 5.0 X lO'^cm and outer radius of nop = 6.72 X lO'^cm, with L = (nop— rhot) = 1-72 X lO'^cm the 
thickness of the spherical shell. Evaluated at mid-depth are the Rayleigh number Ra = (-8 p / dS)(AS / dr)gL^ / pv k, the Taylor number 
Ta = A^'^L^ /u^, the rms Reynolds number Re = Vn-asL/v and fluctuating Reynolds number Re' = v[^^L/u, the magnetic Reynolds 
number Rm = v^^^L/r] and fluctuating magnetic Reynolds number Rm' = v'^^L/r), the Rossby number Ro = uj/IUq and fluctuating 
Rossby number Ro' = uj' /2Qq, and the convective Rossby number Roc = (Ra/TaPr)'/^. Here the fluctuating velocity v' has the 
axisymmetric component removed: v' =v— (v), with angle brackets denoting an average in longitude. The rms velocities in the dynamo 
simulations (and corresponding Reynolds and Rossby numbers) are reduced because the differential rotation is weaker; the fluctuating 
velocities remain comparable. For both simulations, the Prandtl number Pr = j//k is 0.25 and in the dynamo simulation the magnetic 
Prandtl number Pm = u/rj is 0.5. The viscous and magnetic diffusivity, u and rj, are quoted at mid-depth (in units of lO'^ cm-s"'). The 
rotation rate Qq of each reference frame is in multiples of the solar rate ^0 = 2.6 X 10"* rad s"' or 414 nHz. The viscous time scale 
at mid-depth = L^/u is about 3640 days for case D5 and the resistive time scale is about 1820 days, whereas the rotation period is 
5.6 days. For convenient reference, we repeat from Paper 1 the same data for cases D3 and H3 rotating at three times the solar rate. 



Here we focus on another simulation of a solar-type star 
with a faster rotation period of 5 JIq (case D5). As in Pa- 
per I, the simulation builds strong wreaths of toroidal field, 
but unlike those described in Paper I, these wreaths undergo 
quasi-periodic polarity reversals. Three such reversals are 
shown in Figure[Tj7. During a cycle, the global-scale magnetic 
fields wax and wane in strength and can flip their polarity. 
These magnetic activity patterns emerge spontaneously from 
a turbulent convective flow that is significantly more complex 
than in previous laminar, Boussinesq simulations of cyclic dy- 
namos. 

The wreaths of magnetism are highly intricate structures 
with substantial connectivity throughout the convection zone. 
The magnetic wreaths realized in case D5 are shown in field 
line tracings throughout the volume in Figure|2]at a time when 
the magnetic fields are strong. As in Paper I, we find that 
the wreaths are topologically leaky structures, with magnetic 
fields threading in and out of the main flux concentrations near 
the equator Rather than being isolated entities, near the equa- 
tor the two strong wreaths of oppositely directed polarities 
have substantial cross-equatorial connectivity (Fig. This 
cross-equatorial flux appears to play an important role in the 
polarity reversals that are observed. 

Unlike the wreaths of case D3, here magnetism fills the en- 
tire convection zone including the polar regions (Fig. On 
their high-latitude (polar) edges, the wreaths near the equa- 
tor are connected to magnetic structures of weaker amplitude 
and opposite polarity at the polar caps. These polar struc- 
tures are relic wreaths from the previous cycle that propa- 
gate towards the poles during the polarity reversal. In a short 
time after this snapshot, the strong wreaths near the equa- 
tor begin to propagate towards the poles and are replaced by 
new wreaths of opposite polarity (blue in northern and red in 
southern hemisphere). This phenomena is visible in Figure[Tj7 
starting at roughly day 4000, with the reversal completed a 
short time later This is a remarkable example of magnetic 
self-organization by turbulent, rotating, stratified convection 
that bears strongly on the vibrant magnetic activity and cyclic 
variabihty observed in many young, rapidly-rotating stars. 



4. PATTERNS OF CONVECTION IN CASE D5 

The complex patterns of convection of our dynamo and hy- 
drodynamic cases rotating at five times the solar rate are pre- 
sented in Figure |3] Individual convective cells are shown in 
snapshots of radial velocity near the top of the simulated do- 
main in Figures ^ and b for cases D5 and H5 respectively. 
Both cases share strong similarities in their convective pat- 
terns. Owing to the density stratification the convection is 
compressible and the downflows are narrow and fast, while 
the upflows are broader and slower 

Near the equator the prominent cells are aligned north- 
south and propagate in the prograde direction. The strongest 
flows span the entire convection zone; the weaker cells are 
partially truncated by the strong zonal flows of differential 
rotation. Nearer to the poles (above roughly ±45° latitude) 
the patterns are more isotropic. Networks of downflow lanes 
surround upflows and both propagate in a retrograde fash- 
ion. There is less radial shear and most of the convective 
cells span the full convection zone. In the polar regions, the 
radial velocity patterns have a somewhat cuspy appearance, 
with the strongest downflows appearing to favor the westward 
and lower-latitude side of each convective cell. This may be 
a consequence of the strong retrograde differential rotation in 
those regions. 

The convective downflow structures propagate more rapidly 
than the differential rotation that they establish and in which 
they are embedded. In the equatorial band, these structures 
move in a prograde fashion and at high latitudes in a retro- 
grade sense. Individual convective cells typically persist for 
about 10 days, though some have much longer lifetimes. 

The convective structures in case D5 are quite similar to 
those realized in the hydrodynamic case H5 (Figure IS?), 
though there are some noticeable differences, particularly at 
the mid latitudes (around ±30°). In the hydrodynamic case 
there is little radial flow in these regions, as the strong differ- 
ential rotation shears out the convective cells. This region is 
equatorward of the tangent cylinder, an imaginary boundary 
tangent to the base of the convection zone and aligned with the 
rotation axis. For rotating convective shells, it has generally 
been found that the dynamics are different inside and outside 
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Figure 2. Field line tracings of magnetic wreaths in case D5. (a) Snapshot of two wreaths in the equatorial region at day 3880 (time tl), when the magnetic 
fields are strong. Lines trace the magnetic fields and color denotes the amplitude and polarity of the longitudinal field (red, positive; blue, negative); here the 
view is restricted to the equatorial region, spanning roughly ±30° latitude. Magnetic field threads in and out of two oppositely directed wreaths, with substantial 
connectivity across the equator, (b) Side view spanning slightly more than one hemisphere, showing connectivity from equatoiial regions to polar caps. Relic 
wreaths from the previous magnetic cycle are visible at the poles. 



the tangent cylinder, due to differences in connectivity and 
rotational constraint in these two regions (e.g.. .Busse 1970.) . 
The tangent cylinder in our geometry intersects with the stel- 
lar surface at roughly ±42° of latitude. In our compressible 
simulations, we generally find that the convective patterns in 
the equatorial regions are bounded by a conic surface rather 
than the tangent cylinder (Brown et al. 2008). In case H5 the 
strong differential rotation serves to disrupt the convection at 
the mid-latitudes. In contrast, in the dynamo case D5 the dif- 
ferential rotation is substantially weaker in both radial and lat- 
itudinal angular velocity contrasts. As is evident in Figure[3}2, 
the convective cells fill in this region quite completely. 

In our prior hydrodynamic simulations of convection in 
younger suns we r eported on localized nests of convection 
(iBrown et al.ll2008h . with those most prominent at the highest 
rotation rates. Though there is some modulation with longi- 
tude in the equatorial roll amplitudes here, this modulation 
is less extreme in either case D5 or H5 than in our previ- 
ous rapidly rotating simulation s of stellar convection (e.g., 
case G5 in iBrown et al.l 120081) . This difference appears to 
be linked to our background stratification and feedbacks from 
thermal transport near the top of the domain. Here we have 
attempted to reduce the region of influence of the unresolved 
SGS heat flux F„ which carries flux out the top of the domain 
(IBrown et al.l2008i) . This thinner thermal boundary has larger 
gradients and a steeper profile of the background entropy gra- 
dient dS/dr and thus slightly higher Rayleigh numbers and 
radial velocities. In our broader study of rapidly rotating dy- 
namos we have found that strongly localized active nests of 
convection remain possible in dynamo simulations at the most 
rapid rotation rates (fl > 10 fi©). 

The convection establishes a prominent solar-like differen- 
tial rotation, with a fast prograde equator and slow retrograde 
poles. Figure |3]: shows the profile of mean angular velocity 



realized in case D5, averaged in azimuth (longitude) and time 
over a period of roughly 200 days centered on the time of the 
snapshot in FigureO?- The equatorial acceleration is achieved 
by Reynolds stresses and convective transport that redistribute 
angular momentum and entropy, and build prominent gradi- 
ents in latitude of angular velocity and temperature. Radial 
cuts of n indicate that strong radial shear is present through- 
out the lower-latitudes (Figure |3}/). This differential rotation 
is solar-like in the sense that there is a monotonic decrease of 
n from the equator to the pole. Generally, the profiles here of 
n are more cylindrical than those deduced from helioseismol- 
ogy for the Sun, but this is to be expected for more rapidly ro- 
tating stars. This may also be influenced by our omission of a 
tachocline. From studies of solar convection, it is evident that 
the thermal structure of the tachocline with latitude influences 
the differential rotation profiles in the bulk of the convection 
zone. The main effect is to tilt t he contours towa rd a more 
radial alignment (Rempel"2005'; Miesch et al.ll2006l) . 

The differential rotation achieved is stronger in our hydro- 
dynamic case H5 than in our dynamo case D5. This can be 
quantified by measureme nts of the lati t udinal angular velocity 
shear Ai7iat. Here, as in IBrown et al.l (l2008h and Paper I, we 
define Afiiat as the shear near the surface between the equator 
and a high latitude, say ±60°, with 

A^iiat = f^eq " ^^60 , ( 1 ) 

and the radial shear Af^i as the angular velocity shear between 
the surface and bottom of the convection zone near the equator 
with 

Ailj = 51o.97R0 - ^QJ2Rq ■ (2) 

We further define the relative shear as Afliat/rieq- In both 
definitions, we average the measurements of Aft in the north- 
ern and southern hemispheres, as the rotation profile is often 
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Figure 3. Convective structures and mean flows in cases D5 and H5. {a) Ra- 
dial velocity Vr in dynamo case D5 shown in global Mollweide projection at 
O.95R0 with upflows light and downflows dark. Poles are at top and bot- 
tom and the equator is the thick dashed line, while the stellar surface at Rq 
is indicated by the thin sun'ounding line. This snapshot samples day 3880 
(time tl) when the magnetic fields are strong, (b) Companion hydrodynamic 
case H5. Here stronger differential rotation shears out convective structures 
in the mid-latitudes, (c) Profile of mean angular velocity ^(r, 9) for case D5, 
with (d) radial cuts of f2 at selected latitudes, (e) Meridional circulations 
for case D5, with magnitude and sense of circulation indicated by color (red 
counter-clockwise, blue clockwise) and streamlines of mass flux overlaid. 
The profiles shown in (c-e) have been averaged over nearly 2000 days, span- 
ning a full magnetic reversal. 

slightly asymmetric about the equator 

These measurements are quoted for case D5 and H5 in Ta- 
ble |2] The angular velocity shears in case D5 can vary sub- 
stantially in time as the magnetic fields wax and wane in 
strength. As such, here we quote measurements during one 
magnetic cycle with measurements averaged over the entire 
time interval (with label avg and date range indicated) and at 
periods of strongest and weakest differential rotation during 
this cycle {max and min respectively, at indicated times). In 
the hydrodynamic simulation the differential rotation shows 
far smaller time-variations. The global-scale magnetic fields 
realized in the dynamo case D5 feed back on the differential 
rotation and strongly diminish the amplitude of the angular 
velocity shears as compared with the hydrodynamic case H5. 
This results from both a slowing of the equatorial rotation rate 
and an increase in the rotation rate in the polar regions. 



Table 2 

Global Properties of Angular Velocity 



Case 








Epoch 




1.14 


0.71 


0.083 


3500-5500 




0.91 


0.39 


0.067 


3702 




1.43 


0.98 


0.102 


4060 


H5 


2.77 


1.31 


0.192 




D3 


1.18 


0.71 


0.137 


2010-6980 


H3 


2.22 


0.94 


0.246 





Note. — Angular velocity shear in units of /^rad s" , 
with Af^iat measured near the surface (0.97Sq) and 
Afi, measured across the full shell at the equator The 
relative latitudinal shear AQiai/ficq is also shown. For 
the dynamo cases, these measurements are taken over 
the indicated range of days. In oscillating dynamo 
case D5, these measurements are averaged over a long 
epoch (avg), and are also taken at two short inter- 
vals in time when the differential rotation is particu- 
larly strong (max) and when magnetic fields have sup- 
pressed this flow (min). The hydrodynamic case H5 
is averaged for roughly 300 days and shows no sys- 
tematic variation on longer timescales. Measurements 
for cases D3 and H3 rotating at 3 ^0 ai'e quoted from 
Paper I. 

The angular velocity shear realized in case D5 in both lati- 
tude Afiiat and radius Afi,. is remarkably similar in amplitude 
to that realized in our previous dynamo simulation case D3 
(Paper I) even though the basic rotation rate VIq is substan- 
tially faster. This is in striking and marked contrast to our 
hydrodynamic companion cases (H5 and H3) where faster ro- 
tation leads to greater angular velocity contrasts (Table |2]l. 

The differential rotation profiles realized in these rapidly 
rotati ng simulations a re substan tially in t hermal wi nd balance 
(e.g.,lBrun & Toomre 2002; M iesch et al . 2006; Brown et al.l 
120081) . though large departures do arise near the inner and 
outer boundaries where Reynolds stresses and boundary con- 
ditions play a dominant role. Maxwell stresses are signifi- 
cant in the cores of the magnetic wreaths realized in dynamo 
case D5, but these play relatively little role in the global trans- 
port of angular momentum. During a reversal, these Maxwell 
stresses do however transport angular momentum toward the 
poles, giving rise to bands of quickly flowing fluid that share 
some similarities with the torsional oscillations observed dur- 
ing the solar cycle. This behavior will be explored in Sec- 
tion|9l 

The meridional circulation patterns for case D5 are shown 
in Figure [3^. There appear to be three major cells of circu- 
lation in each hemisphere, with several cells in both radius 
and latitude. Strong rotational constraint in regions outside 
the tangent cylinder likely leads to the significantly cylin- 
drical nature of these slow flows. Some flows along the in- 
ner and outer boundaries cross the tangent cylinder and serve 
to weakly couple the polar regions to the equatorial regions. 
These meridional circulations are very similar to the circu- 
lations found in case H5. In both simulations, the merid- 
ional circulations are weaker, slower and more multi-celled 
than those re alized in simulatio n s rota ting at the slower ro- 
tation rates (iBrown et al.l 120081 120101) . The slower flows 
and multi-cellular nature of these meridional circulations may 
have strong imphcations for flux transport dynamos (e.g . , 
iBonanno et al]|2006t iJouve & BrunlllOOTt iJouve et al.ll2010l) . 
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Figure 4. Complex time evolution in case D5 with flips in polarity of magnetic wreaths, (a) Volume-averaged energy densities for kinetic energy of convection 
(CKE), differential rotation (DRKE) and for magnetic energy in fluctuating fields (FME), in mean toroidal fields (TME) and in mean poloidal fields (PME) as 
labeled. Oscillations on roughly 500-1000 day periods are visible in the magnetic energies and in DRKE, though CKE stays nearly constant, (b) Mean toroidal 
field (Bij,) averaged over entire northern and southern hemispheres (labeled) at mid-convection zone (O.85R0). Early in the simulation, opposite polarities 
dominate each hemisphere. Several reversals occur, along with several extreme excursions which do not flip the polarity of the global-scale field. During the 
interval from roughly day 7700 to 10200, the dynamo falls into pecuHar single polarity states, with one polarity dominating both hemispheres. Bracketed interval 
from day 3500 to 5700 spans one full polarity reversal; (c) expanded view of volume-averaged energy densities during this period, and (d) the mean toroidal field 
with same vertical axis scales as in (a,b). Thick labeled tick marks above a,c indicate time samples used in later images. 



5. OSCILLATIONS IN ENERGIES AND CHANGES OF 
POLARITY 

A striking feature of the convective dynamo case D5 is its 
time dependence. This time- varying behavior is readily vis- 
ible as oscillations of the volume-averaged kinetic and mag- 
netic energy densities, as shown in Figure HJ? at a time long 
after the dynamo has saturated and reached equilibration. 

Here the kinetic energy of differential rotation (DRKE) un- 
dergoes factor of five changes on periods of 500-1000 days. 
As DRKE decreases the magnetic energies increase. Moving 
in concert are the mean toroidal (TME) and mean poloidal 
(PME) magnetic energies. The mean poloidal fields appear 
to lag slightly behind the mean toroidal fields as they both 
change in strength. The fluctuating magnetic energy (FME) 
tracks the largest rises in the mean fields but decouples during 
many of the deepest dips. In contrast, the variations in convec- 
tive kinetic energy (CKE) shows little organized behavior in 
time, and appears to change substantially only when the dif- 
ferential rotation is highly suppressed during the period from 
day 7500 to day 8300. The energy contained in the merid- 
ional circulations (MCKE) is weaker and not shown. Though 
it varies somewhat in time, there is not a clear relation to the 
changes in magnetic energies. These energies are defined as 



CKE=ip 



(V0-(V0))^ 



DRKE=-p(y^)2, 



MCKE= -j5({vry + {vey], 



FME=- 



{B,-{B,)f+{Be-{Be)f + 



TME=— 



(3) 
(4) 
(5) 

(6) 
(7) 



PME = 



'+(B„y 



(8) 



where angle brackets will consistently denote an average in 
longitude. 

Magnetic energies in case D5 can rise to be a substantial 
fraction of the kinetic energies. Averaged over the nearly 
16000 days (about 44 years) shown here, the magnetic en- 
ergies are about 17% of the kinetic energies. During individ- 
ual oscillations the magnetic energies can range from a few 
percent of the kinetic energies to levels as high as 50%. The 
kinetic energy is largely in the fluctuating convection and dif- 
ferential rotation, with CKE fairly constant and ranging from 
15-60% of the total kinetic energy as DRKE grows and sub- 
sides, itself contributing between 40 to 85% of the kinetic 
energy. The magnetic energies are largely split between the 
mean toroidal fields and the fluctuating fields, with TME con- 
taining about 35% of the magnetic energy on average, FME 
containing about 61% and PME containing 4%. The roles of 
these energy reservoirs change somewhat through each oscil- 
lation. At any one time, between 10 and 60% of the magnetic 
energy is in TME while FME contains between 30 and 85% 
of the total. Meanwhile, PME can comprise as little as 1 % or 
as much as 10% of the total. Generally, PME is about 12% 
of TME, but because PME lags the changes in TME slightly, 
there are periods of time when PME is almost 40% of TME. 

These results are in contrast to our previous simulations 
of the solar dynamo, where the mean fields contained only a 
small fraction of the magnetic energy (e.g.. lBrun et al.ll2004l 
where TME and PME comprise about 2% of the total). Simu- 
lations of dynamo action in fully-convective M-stars do how- 
ever show high levels of magnetic energy in the mean fields 
([Browning 2008). In those simulations the fluctuating fields 
still contain much of the magnetic energy, but the mean 
toroidal fields possess about 18% of the total throughout most 
of the stellar interior Simulation s of dynamo act ion in the 
convecting cores of A-type stars ( Brun et aTl 12005 ) achieve 
similar results though when fossil fields are included, those 
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Table 3 

Kinetic and Magnetic Energies 



Case 


CKE 


DRKE 


MCKE 


FME 


TME 


PME 


D5'™'- 


1.85 


4.46 


0.006 


0.55 


0.43 


0.048 




1.70 


2.85 


0.005 


0.50 


0.25 


0.062 




1.85 


7.52 


0.007 


0.39 


0.65 


0.042 


H5 


2.27 


34.3 


0.008 









Note. — Volume-averaged energy densities relative to the 
rotating coordinate system. Kinetic energies are shown for 
convection (CKE), differential rotation (DRKE) and merid- 
ional circulations (MCKE). Magnetic energies are shown for 
fluctuating magnetic fields (FME), mean toroidal fields (TME) 
and mean poloidal fields (PME). All energy densities are re- 
ported in units of lO^ergcm"^^. In time-vaiying case D5, these 
energies ai'e averaged over the intervals defined in Table |2] 
while in case H5 an arbitrary 1000 day interval was chosen. 

dynamos can reach states with super-equi partition magnetic 
field strengths and strong mean fields (iFeatherstone et all 
[2009'). In our rapidly rotating suns, the mean fields comprise 
a significant portion of the magnetic energy in the convec- 
tion zone and are as important as the fluctuating fields. In the 
wreath-building dynamo case D3 rotating at three times the 
current solar rate (Paper I) we found that the mean toroidal 
fields contained roughly 43% of the magnetic energy and the 
mean poloidal fields contained about 4%, with the rest in the 
fluctuating fields. 

The global-scale magnetic fields can reverse their polarities 
during some of the oscillations in magnetic energies. This 
is evident in Figure [JJ? showing averages at mid-convection 
zone of the longitudinal magnetic field (B^) over the northern 
and southern hemispheres. Reversals in field polarity occur 
periodically, with typical time scales of roughly 1500 days. 
These reversals appear to happen shortly after peak magnetic 
energies are achieved, but do not occur every time magnetic 
energies undergo a full oscillation. Rather, for each successful 
polarity reversal it appears that several failed reversals occur 
where the magnetic energies drop and the average fields de- 
cline in strength, only to return with the same polarity a few 
hundred days later 

We focus in the following discussion on one such reversal, 
shown in closeup in Figures|4j:,ii and spanning the interval of 
time between days 3500 and 5700. Two reversals occur dur- 
ing this interval, with the global polarities flipping into a new 
state at roughly day 4100 and then changing back again at 
about day 5500. Detailed measurements of kinetic and mag- 
netic energies during this interval are shown in Table |3] 

6. GLOBAL-SCALE MAGNETIC REVERSALS 

6.1. Toroidal Field Reversals 

The nature of the global-scale magnetic fields during the re- 
versal spanning days 3500-5700 are presented in detail in Fig- 
ure |5] Several samples of longitudinal magnetic field are 
shown at mid-convection zone spanning this time period. The 
timing of these samples is indicated in Figure |4] by numeric 
labels and likewise in Figure ^ which shows azimuthally- 
averaged (B^) at this depth in a time-latitude map that spans 
the reversal. The evolution of the mean toroidal field near 
the base of the convection zone is discussed in Section [8] and 
Figure [TOb. 

Before a reversal, the magnetic wreaths of case D5 are very 
similar in appearance to the wreaths realized in our persistent 
wreath-building dynamo case D3 (Paper I). They are domi- 




Figure 5. Evolution of longitudinal field duiing a polarity reversal in 
case D5. (a) Time-latitude plot of (B^) at mid-convection zone, with sam- 
pling times indicated, (h-f) Snapshots of 6^ in MoUweide projection at 
mid-convection zone (O.SSS©) at times indicated. Between reversals the field 
is dominated by the mean component, but during reversals substantial fluctua- 
tions develop, (g — k) Accompanying samples of azimuthally-averaged (B^), 
showing stracture of mean fields with radius and latitude at same instants in 
time. The companion poloidal fields are shown in Figure|2] 
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nated by the azimuthally-averaged component of B^, while 
also showing small-scale variations where convective plumes 
distort the fields (Figure|5f'). At mid-convection zone, typical 
longitudinal field strengths are of order ±13 kG, while peak 
field strengths there can reach ±40 kG. Meanwhile {B^) is 
fairly antisymmetric between the northern and southern hemi- 
spheres (Figure |5^). Shortly before a reversal, the magnetic 
wreaths strengthen in amplitude and become more antisym- 
metric about the equator 

They reach their peak values just before the polarity change 
at roughly day 4000 but then quickly begin to unravel, gain- 
ing significant structure on smaller scales (Figure At 
the same time, prominent magnetic structures detach from 
the higher-latitude edges and begin migrating toward the po- 
lar regions. Meanwhile, {B^) loses its antisymmetry between 
the two hemispheres, with (B^) in one hemisphere typically 
remaining stronger and more concentrated than in the other 
(Figure Is};). The stronger hemisphere (here the northern) re- 
tains its polarity for about 100 days as the fields in the other 
hemisphere (here southern) weaken and reverse in polarity. 
At this point the new wreaths of the next cycle, with opposite 
polarity, are already faintly visible at the equator 

Within another 100 days these new wreaths grow in strength 
and become comparable with the structures they replace, 
which are still visible at higher latitudes (Figures |5]^,/). The 
mean (B^) begins to contribute significantly to the overall 
structure of the new wreaths, and soon the polarity reversal 
is completed. In the interval immediately after the rever- 
sal, small-scale fluctuations still contribute significantly to the 
overall structure of the wreaths, and B^ has a complicated 
structure at mid-convection zone. At this time, the peak mag- 
netic field strengths are somewhat lower, at about ±20 kG. 
As {Bfj,) becomes stronger, the wreaths return to an antisym- 
metric state, with similar amplitudes and structure in both the 
northern and southern hemispheres (Figures |5^,y). They look 
much as they did before the reversal, though now with oppo- 
site polarities. 

The wreaths from the previous cycle appear to move 
through the lower convection zone and toward higher lati- 
tudes. This can be seen variously in the time-latitude map 
at mid-convection zone (Figure |5}3), in the Mollweide snap- 
shots (Figures HJ?-/), as well as in the samples of (Z?^) 
(Figures |5^-fc). This poleward migration appears to be par- 
tially due to a dynamo wave (Section ^ and partially due 
to hoop stresses within the magneti c wreaths and an associ- 
ated poleward-slip instability (e.g.. Spruit & van Ballegooiien| 
IJ982; Moreno-Insertis et al. 1992, and our Section |9]l. Even 
at late times some signatures of the previous wreaths per- 
sist in the polar regions, and are still visible in Figures |5^,y' 
at day 4390. They are much weaker in amplitude than the 
wreaths at the equator, but they persist until the wreaths from 
the next cycle move poleward and replace them. As they ap- 
proach the polar regions, the old wreaths dissipate on both 
large and small scales, for the vortical polar convection shreds 
them apart and ohmic diffusion reconnects them with the relic 
wreaths of the previous cycle. The i7-effect also contributes 
both to the low-latitude generation of the wreaths and their 
high-latitude decay (SectionlSJ. 

Though reversals occur on average once every 1500 days, 
substantial variations can occur on shorter time scales. Here 
at mid-cycle the wreaths become concentrated in smaller lon- 
gitudinal intervals of the equatorial region (as in Figures |5]f;^ 
at day 4780). At other times the mean longitudinal fields be- 
come quite asymmetric, with one hemisphere strong and one 
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Figure 6. The squared magnitude of the antisymmetric, low-latitude toroidal 
field component B^(t) (black lines, left axes) and the unsigned flux across the 
equator 4"f,(0 (blue lines, right axes) as a function of time for two selected 
intervals when the wreaths of opposite polarity occupy each hemisphere: (a) 
t = 1400-6000 (days), and (b) 11,500-16,500 (days). Vertical dotted lines 
indicate polarity reversals, defined as where Ba{t) changes sign. Labeled tick 
marks at top indicate times tl-t5. Generally, the unsigned flux ^[((f) is larger 
when is decreasing. During the interval 6000-11,500 days (not shown) 
the dynamo is in an unusual state in which wreaths in both hemispheres have 
the same magnetic polaiity. 



weak (i.e., during days 4900-5200) before regaining their an- 
tisymmetric nature shortly prior to the next reversal. 

6.2. Connections Across the Equator 

The oppositely-directed wreaths are not isolated entities. 
Rather, they interact through complex magnetic linkages 
across the equator that contribute to their erosion and sub- 
sequent reversal. 

The magnetic linkages between hemispheres can be quanti- 
fied by the net unsigned magnetic flux through the equatorial 
plane: 



2ir 



\Be{r, 6'eq ,cj))\r sin 6*^ drd0, 



(9) 



where 6'eq = 7r/2 and where r\ = rbot and r2 = rtop. In order to 
relate this quantity to polarity reversals of the toroidal wreaths 
we define the antisymmetric component of the mean, low- 
latitude toroidal field as follows 



Ba{t) = V-' 



0„-S Jo 



where 



e„-s Jo 



hi0)B^{r, 0, (j)y sin edrdOdcl) , 

(10) 

r^ sin 9drd9d(j) = — - r^) cos(6l, - S) . 

(11) 
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Figure 7. Evolution of poloidal field during a polaiity reversal in case D5. Shown are the vector potential of the mean poloidal field {A^) in case D5 at times 
tl-t5 (with corresponding days noted). Colors denote amplitude and polarity (red clockwise; blue, counterclockwise); units are G cm. The stellar surface is 
indicated with a solid line, and the poloidal field is extended with a potential field extrapolation to 1.3 Kq (dashed surface). Reversals of poloidal field appear 
to begin in the equatorial region (ci) and spread through the upper convection zone (h,c) before new reversed polarity appears at the equator (d.e). The toroidal 
fields have largely reversed in the equatorial region by time t3 (see Figure|5]/, i). 



Here h{9) is a step function with h{6) = 1 for 9 < 9^, h{9) = -1 
for 9 > 9e, and h{9,) = 0. We set S = 27r/9 (40°) in order to 
focus on the low-latitude wreaths and then plot the squared 
amplitude of this quantity. 

Figure |6] shows ^^'^ ^fl(0 versus time for the two 

intervals in the simulation when cyclic magnetic activity is 
most apparent. The interval spanning days 6000-11,500 is 
not shown as the dynamo had fallen into a peculiar single po- 
larity state and had temporarily stopped showing cyclic be- 
havior; we defer discussion of that interval until Section [TO] 
Labeled tick marks at the top of Figure|6]5! indicate times tl-t5 
shown previously during the reversal in Figure |5](days 3880 
to 4780). 

Examining this interval, we see that the squared ampli- 
tude of antisymmetric mean field attains a peak value near 
time tl shortly before the reversal and then drops to a min- 
imum as the global-scale fields reverse in polarity (reversal 
denoted by vertical dotted line at roughly day 4100). Some of 
this decrease is due to the poleward propagation of the mag- 
netic wreaths. In Figure ^ we see that at mid-convection 
zone, the wreaths of the preceding cycle have largely left the 
region ±40° latitude by day 4100, leaving the wreaths of new 
polarity at the equator then slowly grows in amplitude be- 
fore again attaining a sharp maximum and reversing in sense 
(near day 5500). The unsigned flux through the equator 
lags somewhat behind the mean fields, peaking near time t2 
and dropping to a minimum at approximately day 4200. The 
unsigned flux, measuring cross-equatorial connectivity, does 
not drop to zero during any of the intervals studied here. 

Clearly, the decay and subsequent reversal of antisymmet- 
ric toroidal wreaths is strongly correlated with enhanced mag- 
netic linkages across the equator To quantify this relation- 
ship, we define the rising and decaying phases of a cycle as 
those intervals when the temporal derivative of B^(f) is posi- 
tive and negative respectively. Then we proceed to compute 
the average value of ^l{t) in declining phases relative to rising 
phases. For the two time intervals shown in Figure |6] this ra- 



tio is (a) 2.7 and (b) 1.9. Over the entire extended simulation 
interval (1400-18,252 days) the ratio of unsigned flux in de- 
clining versus rising phases is 1.7. Thus, there is significantly 
more unsigned flux across the equator when the wreaths are 
declining in amplitude or undergoing a reversal. 

6.3. Poloidal Field Reversals 

The evolution of the mean poloidal field can be followed 
by examining how its vector potential (A^) evolves in time, 
where 

(Bpo,) = (B,>f+(Be>0 = Vx (A^0), (12) 

with unit vectors denoted by hats and angle brackets again 
denoting an average in longitude. The vector potential of the 
mean poloidal magnetic field (A^) is shown in snapshots at 
times tl-t5 during the reversal in Figure |7] with color denot- 
ing polarity and poloidal field lines represented by the overly- 
ing contours. A potential field extrapolation has been used to 
follow the poloidal field above the surface out to a distance of 
L3Rq. 

When the magnetic fields are strong (i.e., at time tl, Fig.lTJj) 
the poloidal field is dominated by odd-£ components, with 
significant dipolar and octupolar contributions. The polar re- 
gions have the same polarity (here negative), while the equa- 
torial region has opposite polarity (here positive). The tran- 
sition between poloidal polarities occurs in the cores of the 
magnetic wreaths, where (Z?^) is strong (near ±15° latitude). 

During the course of a reversal, (A^) grows in amplitude 
in the equatorial region. The equatorial polarity expands to 
fill the upper convection zone (times t2-t3, Figs.|7l7,c) while 
in the lower convection zone the dominant polar polarity be- 
gins to disappear. At the same times, toroidal field of the 
new opposite polarity is appearing in this equatorial region 
(Fig.^, d,h,i). Indeed, by time t3 the wreaths of the new cy- 
cle are well established and the wreaths of the previous cycles 
are propagating towards the polar regions. 

The poloidal field of equatorial polarity from the previous 
cycle (positive in the reversal shown in Figure |7} replaces the 



Magnetic Cycles in a Solar-type Star 



11 



poloidal field at the poles, while in the equatorial region op- 
positely directed poloidal field begins to appear at the equator 
(here negative and visible at time t4). This occurs as the new 
wreaths of toroidal field grow in strength and attain a high de- 
gree of axial symmetry with {B^) contributing substantially 
to the structure of the toroidal field. At this point the rever- 
sal is complete, though flux continues to cancel in the lower 
convection zone, leading to the final reversed state shown in 
Figure |7^. 

7. EVOLUTION OF LARGE-SCALE MOMENTS 

The turbulent convection in case D5 gives rise both to 
small-scale, rapidly evolving magnetism, and to slowly evolv- 
ing ordered fields on larger scales. One manifestation of 
the latter is the generation of wreaths of toroidal field, 
whose strength and temporal evolution we have already de- 
scribed. Another is the presence of remarkably strong dipole, 
quadrupole and octupolar components of the poloidal field. 
These low-order moments of the field do not dominate the 
magnetic energy - indeed, near the top of the simulation do- 
main, modes with spherical harmonic degree I = 1-3 typ- 
ically contain no more energy than modes with i up to 
20. But these low-order modes hold particular significance 
both theoretically and observationally: in the Sun, for in- 
stance, the evolution of the global dipole moment is tightly 
linked to the sunspot cycle, with the phasing relationship 
between the two (the sign of the surface dipole reverses 
near solar maximum) serving as an important constraint on 
models of the solar dynamo (e. g., IWang & Sheeleyl 119911: 
iCharbonneau & MacGregodll997b . 

More generally, the low-order modes of the magnetic field 
are important both as diagnostics of global-scale dynamo ac- 
tion and as mediators that regulate the interaction of a star 
with its surroundings. Because higher-order multipoles fall 
off quickly with radius, the dipole mode dominates the mag- 
netic energy distribution at large distances from the stellar 
surface, and so may contribute most to the magnetic "lever 
arm" that determine s how rapi dly a sta r's rotation is braked 
(e.g.. lWeber & Davis 1967; M att & Pudr itz 2008). Low-order 
multipoles also regulate the interaction of pre-main -sequence 
stars with their surrounding accretion disks (e.g., IShu et al] 
11994ft . with broad implications for the truncation of those 
disks and the formation of protoplanetesimals. Such consider- 
ations warrant a careful consideration of the multipolar struc- 
ture and evolution in case D5. 

7.1. The Wandering Dipole 

The lowest-order moment is the dipole, the components of 
which sum vectorally to give a time-evolving quantity with 
both magnitude and direction. Figure |8]illustrates the the evo- 
lution of the dipole vector over the course of one polarity re- 
versal. 

The location of the positive magnetic pole in the lower- 
convection zone is traced in Figure|8p over a period of roughly 
2000 days, displayed on a surface near the base of the convec- 
tion zone. This track was calculated by filtering the radial field 
B, in spectral space (retaining only £ = 1 modes), transform- 
ing back to physical space, and tracking the longitude and 
latitude of the maximum in the resulting field distribution. In 
Figure [Sjj, the latitude of the pole generally stays close to the 
geographic north or south pole, except during reversals of the 
overall polarity. These reversals are the only time during the 
interval studied here that the pole "tips" to latitudes of less 
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Figure 8. Evolution of dipole field in the lower-convection zone during one 
magnetic reversal, with times t2-t5 noted, (a) Location of the maximum in 
the dipole {fi= I) magnetic field (the positive magnetic pole) as measured 
at O.73R0 is traced by the solid black line, with at time t2 and at the 
same depth shown in background, {h) Relative strengths of the axisymmetric 
(m = 0) and non-axisymmetric (m = ± 1 , added in quadrature) terms of the 
dipole. The axisymmetric dipole smoothly evolves from negative to positive 
during the reversal while the non-axisymmetric fields fluctuate continuously 
and show little sign of the global-scale polarity reversal. A second reversal 
occurs after day 5500. 



then about ±45 degrees. The longitudinal position of the pole 
fluctuates more erratically, with no orderly sense of propaga- 
tion. In examining both the longitude and latitude of the l=\ 
maxima and thus including contributions from modes with 
OT = ± 1 as well as m = 0, we are effectively probi ng the evolu- 
tion of the "equatorial dipole" (e.g., Wang & Sheeleyl 1 1 99 Ih 
in addition to the "axial dipole" associated with the m = 
modes. The strength of this equatorial dipole, and the tip- 
ping of the magnetic poles away from the geographic poles, 
are indications of the non-axisymmetric contributions to the 
magnetism. 

A quantitative assessment of non-axisymmetry is provided 
by Figure [8^', which shows the magnitude of the m = com- 
ponent of the dipole and the amplitude of the m = ± 1 dipole 
components added in quadrature during a roughly 2000-day 
interval around this reversal (with second reversal occuring 
after day 5500). As the mean toroidal field in the wreathes re- 
verses polarity (Fig.|5), the axial dipole diminishes in strength 
and changes sign, while the non-axisymmetric pieces fluctu- 
ate erratically and do not seem to sense the global-scale rever- 
sal. 

The erratic evolution of the non-axisymmetric field com- 
ponents, as compared to the relatively smooth evolution of 
the m = field, suggests that the dominant processes gener- 
ating non-axisymmetric fields are somewhat different from 
those responsible for axisymmetric fields. Both the "wander- 
ing of the poles" in longitude (Fig. [Hh) and the fluctuating 
m = ±1 field amplitudes (Fig. [37) would be expected if non- 
axisymmetric field generation were associated mainly with 
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small-scale magnetic features, that add with random phases 
to yield a small (but non-zero) contribution to the dipole. The 
longitude of the pole might then be expected to undergo a ran- 
dom walk on timescales comparable to the convective eddy 
turnover time. This is consistent with what we observe in 
case D5. 

The smooth latitudinal migration of the magnetic pole, on 
the other hand, suggests that the global-scale field reversals 
that occur in case D5 rely on more than just the random ag- 
glomeration of uncorrected small-scale magnetic elements. 
Thus the global-scale axisymmetric m = fields evolve in a 
fashion distinct from the non-axisymmetric components. We 
note that when we try to assess the evolution of the dipole field 
in the fashion shown here at the mid-convection zone or near 
the stellar surface, we find that the dipole there is much more 
variable in time. At these higher radial levels the dipole shows 
several reversals which do not correspond to the global-scale 
reversals of toroidal and poloidal polarity. Instead in the up- 
per convection zone, the higher-order modes of the poloidal 
field more accurately track the reversals occurring throughout 
the convection zone. We consider these in SectionlT2l 



1000 



7.2. Axisymmetric Dipole, Quadrupole and Octupole 

The temporal evolution of the axisymmetric dipole, 
quadrupole, and octupole moments of the poloidal field can 
be assessed by measuring the amplitudes of the m = com- 
ponents of the radial magnetic field Z?,.. These amplitudes are 
defined here as 



V-- 



Q-- 



o= 




B, -COS, e&mOAeAcI) 



B,(3cos2 6l-l)sin6ld6id(/) 



(13) 



(14) 



B,(5cos^6i-3cos6')sin6'd6'd0 (15) 



where the integral so lid angle is at fixed radius (e.g., 
lArfken & Weberfl 19951) . The quantity V is often termed the 
"axial dipole" in solar physics (e.g., Wang & Sheeley 1991). 

The amplitudes T), Q, and O (the £= 1,2,3 components of 
the m = Q field) are shown over an extended interval of about 
13,000 days in Figure |9] Here, the measurements are inte- 
grated over two spherical surfaces, with one near the top of 
the convection zone (Fig.|9]3, at 0.97 Rq) and another near the 
base of the convection zone (Fig. 19^, at Q.13Rq). The oc- 
tupole moment is generally the strongest, particularly when 
the oppositely-directed wreathes are most prominent, but the 
dipole and quadrupole moments are also significant. The 
typical amplitude of these multipoles is comparable to that 
realized in simulations of solar convection with large-scale 
fields (Browning et al. 2006), and stronger than in simula- 
tions of small-scale dynamo action by solar convecti on with 
potential field lower boundary conditions (Brun et al.il2004l) . 
They are also substantially stronger than the surface dipole 
moment observe d on the Sun (e.g., I Wang & SheelevI [19911; 
ISchiiiver & De R osa 2003). Peak pointwise values of the 3D 
radial field Br typically exceed 1 kG. 

Interpreting the evolution of the axisymmetric poloidal field 
is made challenging by the nontrivial depth dependence that 
the various moments exhibit. Near the top of the convection 
zone, the dipole moment V changes sign about twenty times 
over the 13000 days sampled, although some excursions to 
one polarity or the other are short-lived (Fig. |9]5[). These po- 
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Figure 9. Temporal evolution of low-degree moments of the poloidal field. 
Shown are the axisymmetric (m = 0) components of the dipole, quadrupole 
and octupole moments (T>, Q, and O, with £ = 1,2,3 respectively) of the 
radial magnetic field Br at two depths: (a) near the surface at 0.97 Rq and 
(b) near the base of the convection zone at 0.73 Rq. At depth, the octupole 
moment consistently lags the dipole. During the interval roughly spanning 
days 7300-10200 the dynamo enters an unusual symmetric polarity state and 
the quadrupole state becomes strong and does not reverse. After this period, 
normal cyclic behavior resumes. 



larity reversals are separated by fairly irregular intervals of 
300 to 1500 days. The sign of the axisymmetric quadrupole 
component Q flips somewhat more frequently, undergoing 
a few reversals that are not reflected in the dipole moment. 
Generally the dipole and quadrupole track each other well 
in time, with no consistent lag between the two moments. 
The octupole reverses less frequently than either the dipole 
or quadrupole, and when it does reverse it generally lags them 
slightly in time. 

In the lower convection zone (Fig.j^?), the moments evolve 
relatively slowly, with fewer reversals separated by signifi- 
cantly longer intervals of about 1500 days. The octupole mo- 
ment consistently lags a few hundred days behind the dipole. 
The longer evolution times in the lower convection zone could 
plausibly reflect the longer Alfven times deeper down or the 
longer magnetic diffusion time - for our choice of SGS mag- 
netic diffusivity rj, scaling with p"'/^, both timescales vary in 
roughly the same fashion with depth. Convective time scales 
also decrease toward the bottom of the convection zone and 
the ratio of magnetic to kinetic energy increases. Thus, the 
more stable mean fields may reflect less buffeting by convec- 
tive motions. 

The temporal variations in these low-degree moments are 
linked to the reversals of toroidal field, but the relationship 
between the different field components is complex. Near the 
surface, the octupole tracks the deep-seated toroidal fields 
fairly well, but the dipole and quadrupole moments exhibit 
flips which do not correspond to global-scale reversals. In 
the lower convection zone, it appears that the dipolar reversal 
precedes the octupolar reversals. Typically the reversal of the 
dipole occurs shortly after the octupole is near its peak am- 
pUtude (e.g., at about time t4). The octupole then generally 
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reverses when the dipolar moment is near a maximum in am- 
plitude (e.g., about time t5). Most of the reversals in dipole 
and octupole polarity deep in the lower convection zone cor- 
respond to a flip in the sign of the predominant toroidal field 
in each hemisphere, but typically occur only after the wreaths 
have reversed their polarities in the near equatorial regio n (this 
occurs by time t3 in our example reversal). 

When the dipole and octupole are strong, we generally 
see wreaths that are antisymmetric about the equator When 
the quadmpole becomes stronger (e.g., during the interval 
from roughly 7000-10000 days) the dynamo typically ex- 
hibits toroidal fields that are symmetric about the equator 
This indicates that the global-scale toroidal field is likely be- 
ing produced by the stretching of global-scale poloidal field 
by the shear of differential rotation. 

8. A NONLINEAR DYNAMO WAVE 

The latitudinal propagation and orderly reversal of large- 
scale magnetic fields realized in case D5 constitute a striking 
finding of these simulations. Here we examine how field prop- 
agation and reversal might arise, by examining the spatial and 
temporal dependence of some of the mechanisms that act to 
strengthen or weaken fields. 

Spatial propagation and field reversals may be expected 
whenever the processes that amplify or reduce a magnetic 
field - i.e., source terms in the induction equation - are out 
of phase spatially or temporally with the field itself. In mean- 
field solar dynamo models, such spatial propagation can arise 
through the combined operation of the a-effect and Sl-effect 
in the form of a "dynamo wave' ' (e.g.. l s5xl[T976t I Yoshimural 
[19761) . In mean-field solar dynamo models that employ a spa- 
tial separation between poloidal and toroidal source regions, 
flux transport by the meridional circulation, turbulent pump- 
ing, or turbulent diffusion can produce well-defined, non- 
local, phase relationships between the resulting mean fields. 
For example, in many Babcock-Leighton and interface dy- 
namo models the poloidal field at the surface is closely linked 
to the the toroidal field n ear the base of the convection zone 
at a previous tim e (e.g.. ICharbonneau & MacGregorl 119971 
iJouve et al.ll20Tol) . 

In Paper I, we noted that the generation of mean poloidal 
and toroidal field by turbulent fluctuations is not well repre- 
sented by a simple scalar a-effect. This conclusion also ap- 
plies to the simulation reported here. Yet, assessing the spa- 
tial and temporal phasing between mean poloidal and toroidal 
fields and their source terms is nevertheless essential in order 
to understand the physical mechanisms underlying the spatial 
propagation and cyclic reversals exhibited by case D5. 

As in case D3 of Paper 1, the principal source of mean 
toroidal field in case D5 is the fJ-effect, the conversion and 
amplification of mean poloidal field by the mean shear In the 
language of Paper 1, this production term is 



Pms={{B)-V) (v) 
or r aO 



^^^^^iB,)i.,Uco.O{B,)M ^ (16) 



and the principle source of mean poloidal field is the longitu- 
dinal component of the turbulent emf: 



(17) 



The curl of Epi4> contributes to the time derivative of (i?poi) 
(cf. eq.[T2]i. As in previous sections, angular brackets repre- 



sent averages over longitude and primes denote fluctuations 
about the mean, e.g. v' = v- (v). 

The dynamical balances described in detail in Paper I for 
case D3 largely apply also to case D5. The source terms 
( fT6] l and (T% are opposed by ohmic diffusion and by the 
meridional components of the turbulent emf which act to 
fragment and disperse the mean toroidal field. The dis- 
persal of the mean toroidal field by convective motions en- 
compasses the concepts of turbulent diffusion and magnetic 
pumping but it is generally more complex, with intricate 3D 
structure and subtle nonlinear feedbacks. Advection by the 
meridional circulation can also contribute to the time evo- 
lution of the mean toroidal and poloidal fields but we find 
that its role in these simulations is relatively minor. This 
is in stark contrast to Flux-Transport solar dynamo mod- 
els where the advection of toroidal flux by the meridional 
circul a tion regulates cyclic act i vity (e .g., Wang & Sheelev] 
[19911 iDikpati & Charbonneaul fT999l; Kukeretal. 2001 
Nandv & Choudhuri 2001; Dikpati & Oilman 2006; Rempel 
12006; .Jouve & Brua,2007; .Yeates et al..,200 8). In case D5, 
it is repetitive, systematic imbalances between these multiple 
production and dissipation terms which, unlike in case D3, 
gives rise to cyclic behavior 

Figure [Tob exhibits the the mean toroidal field (B^) to- 
gether with (b) the principle toroidal and (c) poloidal source 
terms in the lower convection zone (at 0.73^0) over an in- 
terval of about 2300 days, spanning one full reversal of the 
global-scale field. The poloidal source term Epi is also shown 
in the upper convection zone (Fig. \TU\d . at 0.90/?©) and con- 
tours of Pms are shown over each quantity, to aid in assessing 
phase relationships. 

It is apparent from Figures [TOb.fc that (B^) and Pms trace 
one another closely and largely possess the same sign, high- 
lighting the role of rotational shear in generating and main- 
taining the wreaths. Closer scrutiny reveals more intricate 
phase relationships. The fi-effect precedes the appearance 
of the wreaths, confirming its role as the principle source 
term. Furthermore, it is skewed poleward relative to the 
main toroidal flux concentration. This, coupled with a sub- 
sequent reversal in sign of Pms on the equatorward edge of 
the wreaths, induces poleward propagation. 

Thus the il-effect amplifies the poleward edge of the wreath 
while simultaneously suppressing the equatorward edge. The 
mean poloidal field is generated in the vicinity of the wreaths, 
moving poleward in conjunction with the toroidal flux and 
the toroidal source term Pms Figure [TOb. Thus, the poleward 
propagation may be regarded as a nonlinear manifestation of 
a dynamo wave. It is nonlinear in the sense that the poloidal 
source term Epi is not linearly proportional to the mean field 
and it is a dynamo wave in the sense that spatial propagation 
is induced by means of the relative phasing of poloidal and 
toroidal source terms. Over time, this reinforcement at higher 
latitudes and cancellation at low latitudes contributes to the 
poleward propagation and ultimately to the reversal. Dynamo 
waves like this appea r to have been c aptu red in early solar d y- 
namo simulations bv Gilmanl (119831) and'Gla tzmaieii (Il985h . 

Generally the turbulent emf Epi (Fig.fTOb) tracks the mean 
toroidal field in time and in space, though its sign is sym- 
metric about the equator while (Z?^) is generally antisymmet- 
ric. Closer scrutiny reveals a slight offset toward the poleward 
side of the wreaths where the fi-effect also operates. Thus, in 
contrast to traditional a-fi dynamos, a systematic phase shift 
between {B^) and Epi appears to contribute to the latitudinal 
propagation. 
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Figure 10. Temporal and latitudinal variation of magnetic fields and the effects that generate them, (a) Longitudinally-averaged toroidal field {B^ ) near the base 
of the convection zone (0.73Rq), as image and contour plot in latitude and time, (b) Production of toroidal fields by mean shear Pus- (c) Longitude-averaged 
turbulent electromotive force Epi which generates poloidal vector potential (A^). (d) Same emf £fi shown in the upper convection zone (0.90i?Q). Contours of 
flvis at 0.73Rq are overplotted in all panels for reference (red positive, blue negative). These contours have been smoothed with a 35-day moving time average. 



FigurefTOl/ indicates that the poloidal field in the upper con- 
vection zone is established largely in response to field gener- 
ation in the lower convection zone. In particular, E-pi does not 
begin changing in the upper convection zone until the fields 
in the deep convection zone are already in the process of re- 
versing (times tl-t2). This term continually alters the mean 
poloidal field in the upper convection zone relatively unim- 
peded, in contrast to the lower convection zone where it is 
almost entirely balanced by ohmic diffusion. 

In order to fully characterize the role of the source terms 
in promoting polarity reversals and poleward propagation, we 
must understand them within the context of the mean poloidal 
field evolution shown in Figure |7] The dominant contribu- 
tion to Pms is from the fi-effect, proportional to (Spoi) • V il. 
Since the contours of 57 are primarily cylindrical, Vfi is pre- 
dominantly directed away from the rotation axis (Fig. |3)- It 



is then this component of the mean poloidal field we must 
consider when following the evolution of Pms, namely (Z?s) = 
(-Bpoi) -Vi where s = rsinO is the moment arm. 

The time span between tl through t5 provides a demonstra- 
tion of the processes involved. At time tl, the predominantly 
octupolar configuration of the mean poloidal field (Fig. [TJj) 
coupled with the perfectly conducting lower boundary condi- 
tion yield a positive {B^ (away from the rotation axis) through 
much of the northern hemisphere. This reverses near the equa- 
tor where (Z?s) is directed toward the rotation axis. The as- 
sociated transition between positive and negative (A^) occurs 
across the wreaths. Rotational shear operating on this poloidal 
field structure through the il-effect accounts for the positive 
sign of /"ms through much of the northern hemisphere at time 
tl (FigurefTOt)) as well as the poleward skewness and the low- 
latitude sign reversal. As time proceeds, the neutral surface 



Magnetic Cycles in a Solar-type Star 



15 



(Acj,) = drifts poleward along with the wreaths and gradually 
evolves from a radial orientation to a more horizontal orienta- 
tion by t3 (Fig.|7j:). This is associated with the opening up of 
the poloidal field in the upper convection zone noted in Sec- 
tion |6.3l as low-latitude loops spread poleward. The net result 
is a reversal in (B^) (from positive to negative at mid-latitudes 
in the northern hemisphere), with a corresponding reversal in 
Pms (FigurellOb). 

Rotational shear operating on the cylindrically inward 
poloidal field at the equator at time t3 soon produces toroidal 
field of the opposite sign (negative in the northern hemi- 
sphere). By t4, the low-latitude wreaths dominate the toroidal 
field structure and the reversal is complete (Fig. |5^). Near 
the equator, Epi begins to reverse sign in the lower convec- 
tion zone around time t2 but remains weak while the wreaths 
move poleward. As the new low-latitude toroidal wreaths be- 
come established between t3 and t4, the reversed Epi grows 
in amplitude near the equator, generating mean poloidal field 
of the opposite sense relative to higher latitudes and previous 
times (Fig.lTji). This transition is likely a culmination of the 
evolving 3D (non-axisymmetric) magn etic l inkages across the 
equatorial plane highlighted in Section 16.21 Subsequent evo- 
lution enhances the octupolar structure of the mean poloidal 
field (Fig.|7^), setting the stage for the next reversal. 

Thus the helical nature of the wreaths promotes their pole- 
ward propagation and subsequent polarity reversal. This is 
not to say that the mean fields are simply twisted tori. Rather, 
the maxima and minima of (A^) are displaced relative to (B^) 
such that the local magnetic helicity density of the mean field 
{Atp){B^) changes sign across each wreath. Yet the result- 
ing magnetic topology exhibits regions of oppositely-directed 
(Bs) near the wreath's poleward and equatorward edges. This 
induces poleward propagation by means of a sign reversal in 
Pms as described above. 

A potential alternative interpretation of wreath evolution 
in case D5 is the poleward slip instability whereby axisym- 
metric rings of toroidal flu x drift poleward as a consequence 
of magnetic tension (e.g., [Spruit & van BallegooiienI 119821: 
iMoreno-Insertis et al.l 1199 2V However, if this were occur- 
ring then the poleward propagation of {B^) would be achieved 
by means of a meridional circulation induced by the Lorentz 
force. Instead, we find that the amplitude and structure of this 
advective contribution to the mean magnetic induction are not 
sufficient to account for the observed evolution of (B^), par- 
ticularly with regard to the poleward propagation. 

9. TORSIONAL OSCILLATIONS 

The strong magnetic fields achieved in case D5 couple 
strongly with the global-scale flow of differential rotation. As 
the fields themselves vary in strength, the differential rotation 
responds in turn, becoming stronger as the fields weaken and 
then diminishing as the fields are amplified. These cycles of 
faster and slower differential rotation are visible in the traces 
of DRKE shown previously in Figure |4j3. We revisit here the 
interval explored in close detail in Sections |6]and|8] spanning 
days 3500 to 5700 of the simulation and one full polarity re- 
versal. 

9.1. Bands of Shear 

The angular velocity VL at mid-convection zone is shown 
for this period as a time-latitude map in Figure WTh . Here 
again the timing marks indicate our fiducial times tl-t5. In 
the equatorial regions, the differential rotation remains fast 
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Figure 11. Time-varying differential rotation in case D5. (a) Time-latitude 
map of angular velocity f2 at mid-convection zone (O.SSSq). There are sub- 
stantial temporal variations at both the equator and high latitudes, {h) These 
are accentuated by subtracting the time-averaged profile of f2(r, d) at each lat- 
itude. Visible are poleward propagating speedup structures at high latitudes 
and more uniform modulations near the equator, (c) Corresponding vaiia- 
tions in Af2iat near the surface (upper curve, green) and at mid-convection 
zone (lower, red). 



and prograde, but with some modulation in time. Promi- 
nent structures of speedup are visible propagating toward the 
poles at the high latitudes. These structures are much more 
evident when we subtract the time-averaged profile of for 
this period at each latitude (Figure [12?). These bands ap- 
pear as strong, tilted fast (red) bands extending poleward from 
roughly ±30° latitude. In the northern hemisphere, three such 
speedup bands are launched over this interval. In contrast, in 
the south only two such bands are evident; a third is perhaps 
launched around day 4500, but it does not survive or propa- 
gate. 

Comparing these features with the propagation of magnetic 
fields shown in Figures [Sja andfTOb over the same interval, we 
find that velocity speedup features are well correlated with 
the poleward migration of mean longitudinal magnetic field. 
The velocity features bear some resemblance to the poleward 
branch of torsional oscillations observed in the solar convec- 
tion zone over the course of a solar magnetic activity cycle, 
though on a much shorter time scale here as befits the corre- 
spondingly shorter time between magnetic polarity reversals 
in these dynamo simulations. 

The speedup bands propagate toward the poles relatively 
slowly. In a period of roughly 500 days they travel about 
40° in latitude, and this angular propagation rate of about 
0.08° day"' or about 1.6 x 10~^rads~' is roughly constant at 
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all depths in the convection zone. At mid-convection zone 
(Q.ESRq) this corresponds to a distance of roughly 415 Mm 
and a propagation velocity of about 0.8 Mm day"' or about 
9 ms"'. Near the base of the convection zone (0.73Rq), this 
corresponds to a distance of about 355 Mm and a propagation 
velocity of 0.7 Mm day"' or about 8 ms"' . 

This is considerably slower than the fluctuating latitudinal 
flows associated with the convection which at mid-convection 
zone have peak speeds of ±200 ms"' during this time period. 
The meridional circulations are more difficult to interpret. At 
mid-convection zone they can have instantaneous amplitudes 
of about ±10 ms"', but the fluctuations are large and the time- 
averaged circulations are much slower 1-2 ms"') though 
they are poleward in sense at mid-latitudes. Near the base of 
the convection zone the time-averaged circulations are weakly 
equatorward with amplitudes of about 1 ms"' and thus act to 
resist the poleward propagation of the speedup bands and the 
magnetic wreaths. 

If we do moving 35-day time-averages of the meridional 
flows in the lower convection zone, we see some evidence for 
variations that track the speedup bands. In particular, the flow 
in the core of the wreathes tends to be poleward with an am- 
plitude of up to 1-3 ms"' (up to 5 ms"' without smoothing). 
This poleward flow is likely induced by magnetic tension in 
the toroidal wreathes, analogous to th e mechanism underly- 
ing the poleward shp insta bility (e.g., iMoreno-Insertis et al.l 
119921: iJouve & Brunll2009h . In the absence of rotation and 
convection, the propagation speed associated with the polar 
slip instability is approximately given by the Alfven speed as- 
sociated with the mean toroidal fields. 



va,<p = 



(18) 



In case D5, the propagating toroidal features have an ampli- 
tude of (B^) ^ 3-6 kG in the mid convection zone, implying 
an Alfven speed of about 30-60 m s"' . Toroidal field strengths 
are larger near the base of the convection zone, 15-25 kG, 
yielding larger Alfven speeds, ~ 100-160 ms"', despite the 
higher density (0.19 g cm"^, as opposed to 0.065 g cm"^ in 
the mid convection zone). The much lower poleward flows in 
the wreathes may reflect the inhibiting influence of rotation 
(iMoreno-Insertis et al.ll 19921) and convective pumping. 

In the lower convection zone, the slip-induced poleward 
flow must also operate against the background equatorward 
flow maintained by the convection (Fig. [3^). This gives rise 
to a horizontal convergence of the meridional flow near the 
poleward edge of each wreath and an associated upward flow 
with an amplitude of about 1 ms"', consistent with mass con- 
servation. This is analogous to the recirculating flows around 
shear-generated magnetic flux structures that rise due to mag- 
netic buoyancy (e.g.. lCUne et a l. 2003b a). 

In the idealized polar slip instability, the poleward propa- 
gation of toroidal bands is associated with the formation of 
prograde zonal flows as the fluid within the bands tends to 
conserve its angular momentum. By contrast, the wreathes 
in case D5 are not isolated flux structures, they have a leaky 
topology that allows plasma to escape. This will also influ- 
ence the slip-induced propagation speed but it is likely to ac- 
celerate it rather than decelerate it because prograde zonal 
flows provide gyroscopic stabilization. The wreathes in case 
D5 do induce prograde zonal flow variations as they migrate 
poleward (Fig.fTTTi. as expected from the polar slip instability. 



However, as we shall see in Section |9]2] the zonal flow varia- 
tions in case D5 are largely produced by the Lorentz force, as 
opposed to the Coriolis-induced flows associated with angular 
momentum conservation in closed flux structures. 

The important role of the mean Lorentz force in accelerat- 
ing the torsional oscillations (Section |972] | and the phasing of 
the fi-effect relative to the wreathes (Section |8]l suggests that 
the poleward propagation of the torsional oscillations may be 
mainly attributed to a nonlinear dynamo wave. The phase 
speed of such a wave is governed primarily by the poloidal 
and toroidal source terms in the mean induction equation but 
it may also be linked to the Alfven speed associated with the 
mean latitudinal field. Values of (Be) in the cores of the pole- 
ward propagating bands are roughly 2-3 kG in the mid con- 
vection zone, implying Alfven speeds of about 20-30 ms"'. 
Near the base of the convection zone these values are larger, 
with mean latitudinal field strengths of 6-9 kG and Alfven 
speeds of 40-60 ms"' . 

With the expanded sensitivity of Figure \TTb . we can see 
that the equatorial modulation appears as fast and slow pulses 
which span the latitude range of ±20°. These variations are 
fairly uniform across this equatorial region. The velocity vari- 
ations at the equator do not correspond with the equatorial 
propagating branch of torsional oscillations seen in the Sun 
(iThompson et al.l f2003h . In the Sun, the equatorial branch 
may arise fr om enhanced cooUng in the magn etically active 
regions (e. g.. 'Spruit" 2d03t lRempelll2006l 1200% . 

The temporal variations of the angular velocity contrast in 
latitude Ailiat is shown for this period in Figure [iTl:. At mid- 
convection zone (sampled by red line) the variations in Ailiat 
can be substantial, with large contrasts when the fields are 
strong in the magnetic cycle (prior to tl) and smaller contrasts 
when the fields are in the process of reversing (t2, t3). Near 
the surface (green line) Ai7iat can show even larger variations. 
These near-surface values of Afiiat are reported in Table |2] 
averaged over this entire period (avg) and at points in time 
when the contrast is large (max, at day 3702) and small (min, 
at day 4060). From periods of highest contrast to lowest. Aft 
changes by about 0.5 /irads"', which is a change of roughly 
45% relative to the long running average. This corresponds to 
a total change of about 4% relative to the frame rotation rate 
of 13 ^rads"' (5 VIq). 

9.2. Angular Momentum Transport 

We now turn to discussing the transport of angular momen- 
tum associated with the speedup bands observed in case D5. 
The bands of speedup are clearly associated with the magnetic 
wreaths which are similarly propagating towards the poles, 
and here we examine the physical processes that lead to the 
angular momentum transport that locally speeds up the flow 
of differential rotation. 

Our choice of stress-free and potential-field/perfectly con- 
ducting boundary conditions at the top and the bottom of 
the shell respectively has the advantage that no net external 
torque is applied and angular momentum is conserved. Con- 
vection and magnetism can however redistribute angular mo- 
mentum throughout the shell. In the general MHD case there 
are five processes that serve to transport angular momentum 
(iBrun et al.ll2004l) . These are the Reynolds stresses from fluc- 
tuating flows (which we denote Frs), the meridional circula- 
tions (Fmc), the viscous torque (Fyu), the Maxwell stresses 
from fluctuating fields (Fms), and the magnetic torque from 
the global-scale fields (Fmt)- The latitudinal component of 
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Figure 12. Angular momentum transport in case D5 in tlie northern liemispliere of tlie lower convection zone (0.73^0). (a) Angular velocity Q as image and 
contour plot in latitude and time. Shown at the same depth, with contours of Q overlain to guide comparison, are latitudinal angular momentum fluxes arising 
from (b) Reynolds stresses (Frs), (c) meridional circulations (fMc)i i'l) viscous diffusion (Fyo), (e) Maxwell stresses from fluctuating fields (fMs)> ™d (/) 
magnetic torques from axisymmetric fields (Fmt)- Indicated amplitudes of fluxes ai'e all scaled relative to Fmt, with negative (positive) fluxes in the northern 
(southern) hemisphere representing poleward transport. All fluxes (b — f) have been smoothed with a 35-day moving time-average, but the contours of f2, with 
no temporal smoothing, give a sense of the variations. 



the angular momentum flux F is thus: 



with 



= Frs + -Pmc + ^VD + fivis + -Pmt , 

FRS = prsin6i(veV^), 
^MC = pr sin 6* ( ve ) ( ( V0 ) -I- rior sin 6*) , 
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These five latitudinal angular momentum fluxes are displayed 
for case D5 in Figure [12] along with the angular velocity it- 
self (FigurefTSb). These measurements are taken in the lower 
convection zone where the magnetic wreaths are strong (at 
0.73/?q) and as in B allot et ah (2007) we do not average these 
contributions in time but rather follow their time evolution 
with time-latitude plots. Here we focus on an interval span- 
ning 3700-4200 days when the speedup bands are launched 
during the reversal. We show only the northern hemisphere 
here, but dynamics in the southern hemisphere are similar. 
Contours of angular velocity fl are superimposed on each 
panel to guide the eye. 

The transport of angular momentum is modulated and has 
significant polar branches that are strongly associated with 
the speedup bands. As the magnetic wreaths propagate pole- 
ward, they accelerate the local angular velocity. This is clearly 



seen in Figure [T2b . where the latitudinal shear is reduced 
with the retrograde (blue) regions disappearing as the pro- 
grade (red) bands migrate polewards. The acceleration of the 
high latitude region is largely due to the large-scale magnetic 
torque associated with the axisymmetric fields in the mag- 
netic wreaths {Fmt, Figure [T2]f ). This magnetic torque traces 
the propagation of (B^) quite well. The negative (positive) 
sense of this term in the core of the wreaths in the northern 
(southern) hemisphere indicates that Fmt transports angular 
momentum toward the poles, thus accelerating flow in those 
retrograde regions. 

The Maxwell stresses from the fluctuating fields (Fms, Fig- 
ure \T2h ) also help move angular momentum polewards but 
contribute more weakly than the axisymmetric fields. In the 
equatorial region Fms continues to operate even during the re- 
versal when the axisymmetric fields, and hence Fmt, are very 
small. The Reynolds stresses (Frs, Figure [T2t') generally act 
to transport angular momentum towards the equator but have 
amplitudes lower than either of the Maxwell stresses. In this 
simulation angular momentum transport by viscous diffusion 
(Fyu, Figure fl2b) plays very little role. 

The Coriolis forces associated with the meridional circula- 
tions works in concert with the Maxwell torques to spin up 
the flow in the wreaths {Fmc, Figure[T2b). This spinup largely 
arises from the poloidal flows associated with the wreaths 
themselves, rather than from the time-averaged meridional 
circulations, which are equatorward at this depth. Before the 
speedup bands begin propagating towards the poles, Fmc is 
large and opposite in sense to Fmt, acting to transport an- 
gular momentum towards the equator In the cores of the 
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Figure 13. Extended history of varying differential rotation in case D5. 
(a) Variations of Q(r, 8) at mid-convection zone, (h) Mean longitudinal field 
(B^), at same depth. Poleward propagating speedup structures ai'e visi ble 
during each magnetic oscillation. The time samples used in Figures |5]and[l4] 
are indicated. 

wreaths however, Fmc is converging, which acts to spinup 
the fluid within the wreaths. After the speedup bands have 
launched (e.g., at about day 3850 in the northern hemisphere), 
Fmc changes sense and acts to help maintain the prograde ro- 
tation in the speedup bands. As the wreaths propagate to- 
wards the poles, Fmc is similar in amplitude to Fmt and of 
the same sense. This large and changing Fmc arises almost 
entirely from the second term in equation (l2ll ; the first term 
is much smaller in amplitude (roughly 1% of Fmt)- When no 
time-averaging is applied, Fmc exhibits large fluctuations on 
both short and long timescales, with instantaneous amplitudes 
five or more times larger than Fmt- 

9.3. Sampling Many Magnetic Cycles in Case D5 

The linked variations of angular velocity and magnetic 
fields are shown over considerably longer intervals of time 
in Figures Wih.b. Here too we see the equatorial modula- 
tion over many magnetic cycles and the poleward propagating 
speedup bands. Asymmetries between the northern and south- 
em hemisphere are evident at many times in different cycles. 
The latitudinal angular velocity contrasts exhibit large varia- 
tions during each cycle, with Afiiat near the surface ranging 
from 0.6-1.55 /irads"' during the interval shown here. Suc- 
cessive magnetic cycles can have distinctly different angular 
velocity contrasts, and there are additional long-term modula- 
tions that span many magnetic cycles. These effects may be 
very important for stellar observations of differential rotation. 

These angular velocity variations are consistently associ- 
ated with variations in the magnetism. The mean toroidal 
field (Z?0) is shown over the same interval in Figure [T3t'. The 
poleward propagating magnetic features shown previously in 
Figures |5]3 are evident throughout this longer time sampling, 
now appearing as nearly vertical streaks in {B^), owing to the 
compressed time axis. From day 1500 to 7300, four cycles oc- 
cur in which wreaths of opposite polarity are achieved in each 
hemisphere. After this period, the dynamo explores unusual 
single-polarity states. Here either both wreaths have the same 



polarity (t6) or a single dominant wreath is built (t7). After 
day 10,700 the dynamo emerges from this state and returns to 
building two wreaths of opposite polarity which flip in their 
sense an additional three times as the simulation continues. 

10. STRANGE STATES AND WREATHS OF A SINGLE 
POLARITY 

Though generally wreaths of opposite polarity are built 
in each hemisphere, these oscillating dynamos occasionally 
wander into distinctly different states. This occurs for case D5 
around day 7300. Instead of the two nearly antisymmetric 
wreaths of opposite polarity above and below the equator, the 
dynamo enters a state where the polarity in each hemisphere 
is the same, as shown in Figures [T4b.fe at day 8903. Here 
two wreaths of same polarity occupy the two hemispheres and 
persist for an interval of more than 500 days. The positive- 
polarity reaches average amplitudes of 18 kG while the 
weaker negative polarity structures have average amplitudes 
of only about 3 kG. The azimuthally-averaged profiles of (Z?^) 
emphasize that these wreaths span the convection zone and 
have the same polarity everywhere. During this interval of 
time, the mean poloidal field is predominantly quadrupolar 

The dynamo can also achieve states where only a single 
wreath is built in the equatorial regions, as in Figures [14];, li at 
day 9590. Here a single strong wreath of positive polarity fills 
the northern hemisphere, with {B^) reaching a peak amplitude 
of -1-18 kG. This unique structure persists for about 800 days 
before the dynamo flips polarity and builds a strong wreath 
of negative polarity. The predecessor of this new wreath can 
be seen in profiles of {B^) where a much weaker structure of 
negative polarity is visible in the lower convection zone. 

The strong magnetic fields realized in the single wreath 
states react back on the convective flows. This is evident in the 
accompanying snapshot of radial velocities at mid-convection 
zone (Figure Il4fe). In a narrow band spanning 0-20° lat- 
itude and coinciding with the strong tube, the upflows and 
downflows have been virtually erased. Fluctuations in and 
vg are also very small in this region, and the flow is domi- 
nated by the streaming flows of differential rotation. Within 
the wreath the total magnetic energy (ME) at mid-convection 
zone is locally about 10 to 100 times larger than the kinetic 
energy (KE), while outside the wreath KE exceeds ME by 
factors of roughly 10 to 10** at this depth. We see similar re- 
striction of the convective flows whenever the magnetic fields 
become this strong. 

The differential rotation itself (Figure [l4]/ ") is largely unaf- 
fected by the presence of the strong magnetic wreath. There 
is no clear signature of faster flow down the middle of the 
wreath. Likewise, there is little sign of the structure in profiles 
of the thermodynamic variables P, T, 5, or p, with the mean 
profile instead dominated by latitudinal variations consistent 
with thermal wind balance. 

11. CONCLUSIONS 

In this paper we have explored dynamo action in a solar- 
type star rotating five times faster than our Sun currently does. 
We find that strong dynamo action can occur in the bulk of the 
convection zone in this rapidly rotating sun and that the result- 
ing magnetic fields are organized on global scales into large 
wreaths. Generally, these wreaths have opposite polarity in 
each hemisphere. They are not isolated flux surfaces and in- 
stead are intricately linked across the equator and to the polar 
regions. 
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Figure 14. Strange single-polarity states in case D5 as sampled at times t6 
and t7. (a) Snapshot of at mid-convection zone, showing two strong 
wreaths of the same polarity, (b) Instantaneous profile of (B^) at same time. 
(c) Snapshot of at mid-convection zone at a time when a single wreath 
is formed, (d) Weaker negative polarity structures are visible in profile of 
(B^ ) at same instant, (e) Accompanying snapshot of Vr at mid-convection 
zone, showing flows strongly affected by magnetism. (/) The instantaneous 
differential rotation, shown here as profile of f2(r, d), is largely unaffected by 
the strong wreath. 



This dynamo shows rich time variation, undergoing global- 
scale magnetic polarity reversals roughly every 1500 days. 
During a reversal, the magnetic wreaths near the equatorprop- 
agate towards the poles and are replaced by new wreaths with 
opposite polarity. Reversals occur in both the mean toroidal 
field {B^) and the mean poloidal field. The phasing between 
the two global-scale fields changes with depth in the convec- 
tion zone, with the toroidal field changing first in the lower 
convection zone but the poloidal field reversing first near the 
surface. 

We have analyzed the mechanisms by which the global- 
scale fields are maintained and destroyed in this dynamo. 
Generally, we find that the toroidal fields are built primarily 
by the fJ-effect Pms, where differential rotation stretches the 
global-scale poloidal magnetic field into mean toroidal field. 
This generation term generally is well correlated with {B^), 
but during a reversal Pms appears to reverse in sign before the 
mean fields themselves do. The mean poloidal field is built by 
the fluctuating emf /spi- This mean contribution arises from 
correlations in the fluctuating velocities and magnetic fields. 
Generally, as we found in Paper 1, fipi does not appear to be 
well represented by a simple a-effect where E-pi is given by 
a{B^) . Rather, Epi peaks on the poleward edge of the wreaths. 

During a reversal, the wreaths from the previous cycle prop- 
agate towards the poles. This appears to primarily reflect a 
dynamo wave, with a poleward-slip instability probably con- 



tributing to the dynamics as well. The spatial offset between 
the primary poloidal and toroidal generation terms, Epi and 
Pms respectively, likely contributes to the poleward propagat- 
ing dynamo wave behavior which we see here. 

Accompanying the poleward propagating magnetic struc- 
tures are changes in the global-scale diff'erential rotation. 
Bands of prograde rotating fluid move poleward along with 
the magnetic wreaths. These speedup bands appear to be 
accelerated primarily by magnetic torques arising from the 
global-scale axisymmetric magnetic fields, with an additional 
contribution arising from poloidal circulations associated with 
the propagating magnetic wreaths. The angular velocity con- 
trast Afiiat can vary by more than 45% during a reversal and 
these variations may represent changes of several percent rel- 
ative to the rotating reference frame. As such, we may ex- 
pect similar variations in observations of surface differential 
rotation on magnetically active stars during the course of a 
magnetic cycle. The bands of speedup bear some resem- 
blance to the poleward branch of the torsional oscillations 
see n in observations of the Sun over the course of a solar cy- 
cle ([Thompson et al.ll2003h . In these simulations, we do not 
see any evidence of the observed equatorial branch. In the 
Sun, the equatorial branch may arise from enhanced coolirig 
in the magnetically active regions (e.g.. Spruit 2003; Rempe| 
120061 12007) . but those effects lie beyond the spatial resolution 
or physical effects included in our simulations of global-scale 
convection and dynamo action. 

These global-scale polarity reversals are not special to 
case D5. Indeed, we have explored a broader class of oscillat- 
ing dynamo solutions, which will be detailed in forthcoming 
papers. Some of these solutions are realized by taking our 
more slowly rotating case D3 (Paper I) to higher levels of tur- 
bulence by reducing the eddy diffusivities, while others are 
achieved in simulations spinning at even higher rotation rates 
than case D5. Large oscillations in the magnetic fields and 
global-scale reversals of polarity appear to be common fea- 
tures in the parameter space we have explored, likely arising 
when the ohmic diffusivity on large scales is small enough to 
allow the dynamo waves to run polewards in these wreath- 
building dynamos. We find such global-scale oscillations and 
polarity reversals fascinating, since these are self-consistent 
3D stellar dynamo simulations which achieve such temporally 
organized behavior in the bulk of the convection zone, with- 
out appealing to a stable tachocline of shear at the base of the 
convection zone as the organizational seat of the dynamo. 
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